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Abstract. Superconductivity in the iron pnictides and chalcogenides arises at the 
border of antiferromagnetism, which raises the question of the role of quantum 
criticality. In this topical review, wc describe the theoretical work that led to the 
prediction for a magnetic quantum critical point arising out of a competition between 
electronic localization and itinerancy, and the proposal for accessing it by using 
isoelectronic P substitution for As in the undoped iron pnictides. We go on to compile 
the emerging experimental evidence in support of the existence of such a quantum 
critical point in isoelectronically-tuned iron pnictides. We close by discussing the 
implications of these results for the physics of the iron pnictides and chalcogenides. 
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1. Introduction 

In recent years, stimulated by the discovery of new materials with unusual properties, the 
physics of "quantum criticality" has become a frontier issue in condensed matter physics. 
In many materials, several ground states are often in competition and it is usually 
possible to tune from one ground state to another by adjusting some control parameter. 
Such a transition from one zero temperature phase to another is called a quantum 
phase transition; when the transition is continuous, it defines a quantum critical point 
[1]. Typical control parameters include magnetic field, pressure, electron density, and 
chemical composition. In contrast to ordinary thermally driven phase transitions, where 
the control parameter is temperature, dynamical fluctuations of order parameters are 
on equal footing with the usual spatial fluctuations. The quantum critical scaling that 
describes the neighborhood of quantum critical points has been widely discussed starting 
from the very early work of John Hertz |2]. A recent review (in the context of Fermi- 
liquid instabilities) [3] summarizes important aspects of quantum criticality. Just as 
spatial fluctuations are characterized by a correlation length ^ that diverges at a phase 
transition governed by a control parameter x, 



dynamical (temporal) fluctuations are characterized by a diverging correlation time 



Here, z is the dynamic exponent of the quantum phase transition. 

Examples of well-studied quantum phase transitions are the zero temperature 
transitions in quantum Hall effect systems [1], in heavy fermion compounds 0, in the 
cuprate high-temperature superconductors [6], and, in disordered metals, the metal- 
insulator transition [7]. Theoretical work has advanced the notion that Landau's 
paradigm of order-parameter fluctuations, which works so well for thermally-driven 
phase transitions, can fail for quantum critical point [U |9]. Consequently, a unified 
theory of quantum criticality does not yet exist, and many fundamental questions 
are open, for example: What collective excitations are generally important at and 
near a quantum critical point? Are there universality classes? What are the scaling 
laws and their exponents? Under the circumstances, the investigation of different 
classes of quantum critical phenomena in new materials can be an avenue to further 
understanding. 

The discovery [TOl |TTj of a number of iron-based compounds that exhibit 
superconductivity at relatively high temperatures (compared to "conventional" 
superconductors) has generated a large body of experimental and theoretical research 
[121 [13] and because of the presence of competing phases in these compounds we have 
a new window into quantum critical phenomena. 

The phase diagrams of the iron pnictides rather universally exhibit an 
antiferromagnetic ground state for the undoped (parent) compounds, e.g. LaFeAsO 
(1111), BaFe2As2 (122). Upon substitution of F for O in the 1111 case and K for Ba or 
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Co for Fe in the 122 case, the ordered moment gradually disappears as a critical doping 
level is approached. However superconductivity intervenes in this doping region and 
magnetic quantum criticality is not so evident. In figure [l| the behavior of Co doped 
122 is shown [T5] . which illustrates this point. Therefore it would be advantageous 
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Figure 1. Phase diagram of Ba(Fei_2.Coa;)2As2. Adapted from Wang et al [H] 
and Chu et al [15]; the antiferromagnetic (AFM) state is also structuraUy distorted. 
In this case, the quantum critical point at i « 0.2 is obscured by the onset of 
superconductivity with a maximum transition temperature of ~ 25 K. The dashed 
line is the expected continuation of the magnetic transition line in the absence of 
superconductivity. 

to seek an avenue for the exploration of the magnetic quantum critical point that is not 
obscured by the infiuence of other phases. 

The influence of electron-electron interactions, i.e. the consequences of correlation 
effects, plays a central role in our discussion and we begin with remarks that review that 
issue. Then we discuss some consequences of proximity to a Mott transition. This is 
followed by considerations of a magnetic quantum criticality that is tuned by electron- 
electron interaction strength rather than by the usual parameters mentioned above. We 
give examples in the iron-based compounds and conclude with a summary and outlook. 

2. Iron pnictides as bad metals 

The materials under discussion in this review are paramagnetic metals at temperatures 
greater than about 200 K and usually become antiferromagnets at lower T. When 
doped away from their stoichiometric compositions the magnetism disappears and 
superconductivity often arises. A recent review ^12j summarizes the physical properties 
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of the various compounds. Soon after the discovery of superconductivity in the 
iron arsenide compound, we pointed out p!6] that electron-electron interaction has 
a significant effect on these properties. Subsequently, we emphasized and reviewed 
comprehensively [17] the evidence for the importance of correlation effects. In this 
section, we briefly summarize those considerations. 

What are diagnostics to evaluate the strength of electron-electron correlation in 
a metal? One is the magnitude of the room temperature resistivity. Electron-phonon 
scattering rarely leads to a room temperature mean free path i that is of order the lattice 
spacing. Therefore, a small value of kpi can be an indicator of strong electron-electron 
interaction. In the pnictides, there are five bands that cross the Fermi surface. For a 
rough estimate, let us assume that the corresponding Fermi pockets all have the same 
value of kpi. The in-plane conductivity (TB] is the sum of contributions of each pocket, 
say p of them: 

cr = p{e'^/h)kFi/c, 

where c is the c-axis lattice spacing. A typical value for the single-crystal in-plane 
conductivity at room temperature is 2500 (Ohm-cm)~^ and c is about 6.5 A for 122 
iron arsenides and about 8.5 A for 1111 iron arsenides. We get kpi ^ 5/p, which gives 
kpi of order unity for each pocket. Hence, the pnictides may properly be described as 
"bad metals" as a consequence of strong electron-electron interactions. 

Another indication of the importance of correlation comes from a diminished 
strength of the Drude peak of the optical conductivity. A particularly comprehensive 
discussion of this maybe found in a paper of Qazilbash et al [IH] (see also, reference [121 
|20]). in which the pnictides are compared to well-known strongly and moderately 
correlated materials such as the cuprate superconductors and V2O3. The reduction of 
the Drude weight provides a measure of the distance of the bad-metal system from the 
Mott transition point. By this measure, BaFe2As2, for instance, which has only about 
30% of the electronic excitations residing in the coherent part of the optical spectrum 
[18] . has interactions that are about 70% of the Mott-transition value. 

The reduction of Drude weight is accompanied by unusual features in the 
temperature-induced transfer of spectral weight. It is observed [21], [221 ESj that as 
the temperature is lowered, there is a transfer of optical spectral weight from low to 
high (on the order of eV) energies. This is expected as a precursor to the formation of 
Hubbard bands [24j and Mott-Hubbard localization. 

3. Mott transition 

If the electron-electron interaction is sufficiently strong, it is expected that the ground 
state of the system will be a Mott insulator, characterized by upper and lower Hubbard 
bands separated by an energy gap - the so-called Mott gap, of order the interaction 
strength minus the bandwidth. There are many examples, notably the undoped cuprate 
superconductors and vanadium sesquioxide. The ratio of the interaction strength to the 
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bandwidth, i.e. the ratio of potential to kinetic energy, as usual, measures the proximity 
to the Mott localization transition. On the localized side of the transition, the low 
lying excitations are of the degrees of freedom associated with localized spins, whose 
spectral weight comes from the lower Hubbard band. We have argued jJH [17] that the 
iron pnictides are on the metallic side of the Mott transition, but very close to it. In 
this picture, the electron spectral function has a three peak structure consisting of the 
precursors of the Hubbard bands away from the Fermi level Ep and a quasiparticle peak 
at Ep that is responsible for the metallic behavior, as depicted in figure |2j It is then 
natural to ascribe [16] the magnetism of the undoped metallic iron-based compounds to 
the behavior of the local moments associated with the developing lower Hubbard band. 
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Figure 2. Single-electron spectral function as the sum of coherent and 
incoherent parts. The single-electron density of states (DOS) is plotted against 
energy (E); Ep is the Fermi energy. Each peak may contain additional structure due 
to the multi-orbital nature of the iron pnictides. 

The situation depicted in figure [2] is a simplified sketch of what has been determined 
from electron structure studies by a Rutgers University group that has taken correlation 
effects into account systematically |2S]- In that work, in analyses of the optical 
conductivity, the authors find a strongly reduced Drude peak and substantial incoherent 
spectral weight smeared over energies of the order of an electron volt (see also, 
reference |26]). The absence of well-defined features associated with incipient localized 
3d states is due to their crystal field splittings as well as broad features from the As 4p 
states. Their analysis also indicates that although the system is metallic, the sublattice 
magnetization with reduced moment resides on the iron atoms. 

The standard theoretical description for Mott transition occurs in one-band models 
with an odd number of electrons per site. The parent iron pnictides and chalcogenides 
contain multiple orbitals that are partially occupied by an even number of (six) 3d 
electrons per Fe site, and any insulating state could in principle simply be a band 
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insulator. Multi-orbital models for the iron pnictides have been studied in Ref. [27] 
with the conclusion that a Mott transition robustly exists as a function of U/ 1. 

The description of the properties near the Mott transition suggests that by tuning 
the ratio of potential to kinetic energy ("PE/KE") one might be able to access two 
quantum phase transitions. One is the Mott metal-insulator transition itself; the other 
is a magnetic transition from antiferromagnet to paramagnet on the metallic side as one 
moves away from the Mott state. In this way, we interpret the observed magnetism in the 
undoped arsenides as due to quasilocalized spins belonging to the precursor of the lower 
Hubbard band, as shown in figure [2j In this paper we begin with the magnetic transition 
and discuss the Mott transition near the end in connection with recent experiments on 
iron oxychalcogenides. 

The principal method of tuning PE/KE that we consider here is by isoelectronic 
chemical substitution that affects the lattice parameters. The advantage over charge 
doping {e.g. F for O, K for Ba, Co for Fe) is that neutral impurities create much less 
disorder than charged ones. So much so that deHaas van Alphen experiments succeed 
over a substantial range of P doping for As (see section 5.1) and NMR lines remain 
narrow [28J. Specifically, upon substitution of smaller phosphorous for larger arsenic 
in an iron arsenide, the lattice constant decreases; therefore the electron bandwidth, 
the KE in the tuning ratio, increases. At the same time, the intra-atomic interactions 
responsible for correlation effects remain essentially constant. Thus, one may pass from 
the antiferromagnetic state of the arsenide to the putative paramagnetic state of the 
phosphide via a quantum critical point at a critical concentration of phosphorous. This 
was the prediction made in |30] and subsequently confirmed experimentally [31]. The 
evidence is reviewed in section 5. 

4. Theory of the magnetic quantum critical point 

The antiferromagnetism in a typical iron pnictide is characterized by Q = (0,7?) or 
(tt, 0) order on a square lattice [31]. That is, in one direction there is antiferromagnetic 
order and in the orthogonal direction, ferromagnetic order. For the present, we ignore 
the orthorhombic distortion that accompanies the magnetism. Consistent with the 
discussion of the previous section, in which the magnetism arises from the localized 
spins in the precursor Hubbard band, the exchange coupling naturally involves first 
(Ji) and second neighbor (J2) spins with J2 > Ji/2 [16]. The magnetic part of the 
Hamiltonian is then [32] 



This model, with the prescribed ratio between J2 and Ji, naturally generates the 
collinear order [331 El] , and was proposed to underlie the magnetic order in the parent 
iron pnictides [16] . 

It is worth underscoring the preconditions that make this spin Hamiltonian the 
starting point of the description for the spin properties of the parent iron pnictides 
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[161 [30] . The considerations in the previous two sections imply that the spin spectral 
weight will primarily come from electronic states away from the Fermi energy. We use 
a parameter w to measure the coherent part of the electronic spectral weight near the 
Fermi energy (c/. figure [2]). w is relatively small when the electronic correlation is 
close to the threshold value for the Mott transition. It is also small because the Fermi 
surfaces of the parent iron pnictides (and chalcogenides) comprise small electron and 
hole pockets. The weight w can therefore be used as the perturbation parameter, in 
place of U/t that is order unity. The spin Hamiltonian ^ arises in the zero-th order 
of w. The effect of higher order terms on the spin properties is treated below; cf. eq 
nation 




sublattices 
■ A 

♦ B 

• C 
A D 



Figure 3. Classical ground state of the Ji — J2 model. Two interpenetrating 
Neel square lattices with staggered magnetizations m and m', where </) is the angle 
between them. 

As discussed in |i34j, when an antiferromagnetic J2 exceeds Ji/2, the classical ground 
state consists of two interpenetrating Neel lattices (the dotted squares in figure [s]) with 
independent staggered magnetizations (Neel vectors) m and m'. Thus, the order is 
characterized by four sublattices imbedded in a square lattice, as shown in figure|3j Here, 
sublattices A and B are antiferromagnetically coupled with staggered magnetization m 
as are C and D with m'. Within the Ji — J2 model, the energy is independent of the 
angle (p between m and m'. The ultimate coUinear order in the quantum system is 
achieved by the phenomenon called "order from disorder" [35]. Here, the "disorder" 
refers to thermal or quantum fluctuations [33] that break the degeneracy in and 
favor the collinear arrangement = or vr. Thus m ■ m' = ±1 (or equivalent ly a, 
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where m' = am [36]) becomes an Ising order parameter and we shall include it in the 
Ginzburg-Landau theory that we develop below. The identification of an Ising transition 
in the frustrated Ji — J2 model was first made in |3l] and that formulation was applied 
to the pnictides in [361 EZ] • 

The Ising order breaks the C4 symmetry and naturally leads to a tetragonal-to- 
orthorhombic distortion in the presence of a coupling between the Ising order parameter 
m • m' and a structural degree of freedom. The structural distortion, perhaps in 
combination with orbital ordering among the Fe dxz and dyz states [38l 139] . can give 
rise to a planar anisotropy in the value of Ji, i.e. Jia 7^ Ju [38], as has been suggested 
by a spin-wave analysis of the inelastic neutron-scattering spectrum of CaFe2As2 in its 
magnetically ordered phase 

Equation ^ applies at the Mott transition point, where w = 0. The effect of a 
non-zero w can be treated order-by-order in terms of a Ginzburg-Landau action. As we 
discussed in the previous section, the quasiparticle weight w at the Fermi energy seen in 
figure |2] can be tuned by the strength of electronic correlation. We propose that when 
w is increased, a quantum phase transition out of the coUinear antiferromagnetic state 
occurs. We consider the lattice spacing as a tuning parameter but the considerations 
here just depend on the spectral weight w of the quasiparticle peak at the Fermi energy: 
the smaller the lattice parameter, the larger is the KE (or, equivalently, the smaller 
is U/t) and the larger is w and we shall use it as a tuning parameter. With all these 
considerations, the quantum Ginzburg-Landau action at T = has the form [TTf [30] 

5(m,m') = J d{q} J dM[S2(q, u;) + S4({q}, M) + . . .], 

S'2(q, u) = [ro + wAq + cq^ + u'^ + 'j\uj\] [|m(q, cu) |^ + |m'(q, uj)\^] 

+v{ql - gj)[m(q,a;) • m'(-q, -w)], 
S'4({q}, {to}) = "udml^ + |ni'|^) + n|m ■ m'p + -u'lml^ |m'| ^. (3) 

The q-dependent terms in ^2 appear naturally from a gradient expansion of the energy 
represented in (|2]). The notation {q}, {a;} denotes the set of four momenta and four 
frequencies that enter 5*4. The momentum and frequency integrals are understood to 
each contain a delta function that fixes the sum of the momenta and the sum of the 
frequencies at zero. Several important comments on ^ are in order: 

We take ro < 0, so that when w = 0, the system is antiferromagnetically ordered. 
At a critical value Wc, tq+WcAq = 0, which defines the quantum critical point separating 
the magnetic state from a non-magnetic one. A sketch of the resulting phase diagram is 
shown in figure |4} The dashed red crossover lines border the region of quantum critical 
fiuctuations and are defined by h/ksT ~ r, or T oc |^y — WcT^, see (lb). 

When w is non-zero, the local moments are coupled to the itinerant quasiparticles; 
this leads to the damping term ■jIujI in S2, where for low energies, 7 is 0(w°) ^30]. This 
is valid as long as |ci;| is less than the quasiparticle bandwidth wW. Here W is the full 
bandwidth at w = 1. Therefore, the critical theory has dynamical exponent z = 2 [IT] , 
as is typical for an antiferromagnetic quantum phase transition. 
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Figure 4. Pnictide phase diagram near a magnetic quantum critical point. 

The red dot denotes the quantum critical point determined by the critical tuning 
parameter Wc- The blue line is the line of the thermally-driven antiferromagnetic 
transition and the dashed line is a structural transition. At non-zero temperatures, the 
two transitions can occur either coincidentally, or be separated. At zero temperature, 
they coincide. 

As mentioned earlier, the non-zero temperature magnetic transition in the parent 
iron pnictides is accompanied or preceded by a tetragonal to orthorhombic structural 
distortion in which the square array of iron atoms becomes rectangular, where the longer 
a-axis is the direction of the eventual antiferromagnetic order. This can naturally be 
accommodated in the theory by including in (jsj) a coupling of the Ising order parameter 
m ■ m' to a structural degree of freedom. 

The sign of the quartic coupling m in ([3]) is very significant. It can be derived 
from the above Ji — J2 model |34j; it is negative, reflecting the physics of order from 
disorder inherent in the Ji — J2 model [M] and therefore leads to the Ising transition 
fixing m ■ m' = ±1 as discussed above. It also gives rise to the coUinear arrangement 
in the ordered state, in which the moments are antiferromagnetically ordered along the 
a-axis and ferromagnetically ordered along the 6-axis of the planar Fe-Fe lattice. 

In the traditional approach to quantum criticality [2l |12] it is recognized that 
deff = d + z is the effective dimensionality of critical fluctuations, where d is the spatial 
dimension. In our d = 2, z = 2 case, the deff is 4; therefore we are at the upper 
critical dimension at which there are logarithmic corrections to simple Gaussian critical 
behavior and some predictions can be made from scaling arguments [3]. Experimentally, 
the quantum critical behavior can be explored either by tuning w through Wc at very 
low temperature or by lowering the temperature with w fixed at Wc- For example, the 
specific heat varies as —T\n{w — Wc) and — TlnT respectively for the two cases. The 
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resistivity is determined by the interplay of two scattering mechanisms: spin fluctuations 
and impurities. In the dirty limit, for w = Wc, the resistivity varies linearly with T ^3] , 
p — Po ocT, rather than quadratically as in a Fermi liquid. In this d = z = 2 case, the 
uniform susceptibility is expected to vary as x — Xo oc (— TlnT) 

5. Accessing the quantum critical point by isoelectronic substitution 

We mentioned at the end of Section 3 that substitution of P for As decreases the lattice 
constant, thereby reducing the PE/KE ratio and increasing w. Table [l| gives some 
representative numbers that illustrate the change in lattice parameter. The references 
are for the La and Ce arsenides {x = 1), for the La and Ce phosphides, [TF] for 
Ba, and |18] for the Eu compound. 



Table 1. Lattice parameters for isoelectronic substitution. 



Compound 




X 


c(l) 


a(A) 


c/a 




Reference 


LaFe(As2.Pi_^ 


.)0 


1 


8.754 


4.038 


2.168 


142.7 


m 









8.508 


3.957 


2.150 


133.2 


m 


CeFe(As:rPi_ 


.)0 


1 


8.656 


4.000 


2.164 


138.5 











8.328 


3.919 


2.125 


127.9 


m 


BaFe2(As^Pi_ 


-x)2 


1 


13.001 


3.962 


3.281 


204.1 











12.422 


3.844 


3.232 


183.5 


[47^ 


EuFe2(As^Pi. 


-2)2 


1 


12.136 


3.910 


3.104 


185.8 


mi 






0.7 


11.831 


3.889 


3.042 


178.9 





It can be seen that isoelectronic P substitution for As results in a substantial 
reduction in unit cell volume {V), which implies an increase in KE with little change in 
PE, which is primarily local. Thus as x decreases, the tuning parameter w increases. 
Since the As parent compounds exhibit antiferromagnetic order (small w side of figure 
3) and the P compounds do not, we can place the two end materials on the two 
sides of Wc along the tuning-parameter axis in the theoretical phase diagram, figure 
[ij Consequently, one expects at some value of x = Xc (corresponding to Wc) to access 
a magnetic quantum critical point. Then the measurement of various properties at P 
concentration Xc will give information on the quantum critical behavior, e.g. the critical 
exponents. Thus, the method of isoelectronic substitution offers a new route to quantum 
criticality, as distinct from charge carrier doping or application of external pressure. At 
the same time, this approach is not necessarily free of the interference of other phases 
as discussed above in connection with figure [T} 

5.1. 122 compounds 

These are the iron pnictides of composition QFe2Pn2, where Q is usually Ba, but can also 
be a different valence 2 element e.g. Ca, Sr, and even Eu. To achieve superconductivity. 
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the 122s are either electron doped by partial substitution of Fe by another transition 
metal ion, or hole doped by partial substitution of Q by an alkali metal ion. However, 
P doping for As in e.g. BaFe2As2 also causes the suppression of magnetism and the 
appearance of superconductivity. In spite of the interference of superconductivity 
at the critical P concentration Xc for the disappearance of antiferromagnetism, some 
information has been obtained about transport and thermodynamic properties at x = Xc 
in (Ba or Eu)Fe2(Asi_^P^.)2. 

In EuFe2(Asi_2,'Pa;)2 with x = 0.3 |48], no trace of the AFM order remains and not 
only is there an apparent superconducting transition at T = 26 K, but the Eu moments 
appear to order ferromagnetically at T = 20 K. Nevertheless, at higher temperatures, 
the resistivity is linear, p = po + aT, in temperature from about 90 K up to 300 K. 
However, the specific heat does not show signs of the T In T behavior expected for the 
d = 2, z = 2 described above in section 4. While the appearance of the two phases 
that are induced by P-doping is interesting in itself, it obscures the magnetic quantum 
criticality associated with the suppression of the antiferromagnetic state. 

A group at Zhejiang University has studied the BaFe2(Asi_a;Px)2 series [19]. The 
experimentally determined phase diagram is very similar to that shown in figure [TJ The 
principal difference is that the scale of doping should be rescaled by a factor 1.6. The 
fits to a resistivity of the form p = Po + clT'^ give a near 2 (as in a Fermi liquid) for 
concentrations x far from Xc but near the critical concentration Xc = 0.32, a T-linear 
resistivity was found up to room temperature: a ~ 1 for 0.3 < x < 0.45. 

The T-linear behavior of the resistivity in the critical concentration region of 
BaFe2(Asi_3;Pa;)2 was subsequently confirmed by a group from Kyoto University [50] . 
In fact, as emphasized in [SD], the region of the phase diagram in which the resistivity is 
linear in T bears a striking resemblance to the quantum critical region shown in figure 
|4j in spite of the presence of superconductivity at the critical concentration. 

An examination of the temperature dependence of the amplitude of oscillations in 
the de Haas- van Alphen (dHvA) effect can give information on the effective mass(es) 
associated with the extremal portions of the Fermi surface(s) [SI]- In the iron pnictides, 
both electron (at the Brillouin zone edges) and hole (at the zone center) Fermi 
surfaces (pockets) of the essentially two-dimensional electronic structure have been 
identified both from band structure calculations and from experiment. H. Shishido 
et al, a multi-institutional collaboration, [53] have performed dHvA experiments on the 
BaFe2(Asi_2;Pa;)2 series. By extracting the temperature dependent amplitude factor 
{X/ sinhX, X oc m*T / B) that appears in the Lifshitz-Kosevitch formula [5l] for the 
oscillations, they have determined the x dependence of the effective mass on two of the 
pockets. They argue these must be electron pockets at the Brillouin zone edge because 
the hole scattering rate is probably too large to allow oscillations on the zone-center 
hole pockets to be detectable ^52j- At low temperatures (T < 4 K), m* appears to 
diverge as the concentration approaches Xc ~ 0.33 in a manner that is not inconsistent 
with a quantum critical point. It is important to recognize that for essentially all known 
superconductors (in zero magnetic field) the specific heat shows a mean-field-like jump. 
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which indicates that the effect of classical critical fluctuations on the mass is minimal. 
Thus the observed behavior here may be ascribed to the presence of the magnetic 
quantum critical point beneath the superconducting dome. As mentioned at the end of 
section 4, the specific heat, hence m*, should vary as — \n{w — Wc). So we may expect 
the experimental effective mass at low temperature to have a form a — bln{x — Xc)- 
Figure |5] shows this fit for one of the pockets for the few points above Xc- 




Figure 5. Effective mass of BaFe2(Asi_2.P2.)2. The experimental data (black 
dots) is taken from figure 3b of reference 53 . A fit ni* /m — 0.88 — 0.89 ln(a; — 0.33) 
is shown by the solid line. 

The group at Kyoto University has also carried out nuclear magnetic resonance 
(NMR) experiments on the BaFe2(Asi_^.P2.)2 series [281 129]. For the '^^P nucleus, they 
examined the behavior of nuclear spin-lattice relaxation rate 1/Ti and Knight shift K as 
a function of x and T. K depends on the uniform (q = 0) electron susceptibility xIq); 
while 1/Ti is determined by the imaginary part of the local (q-averaged) susceptibility. 
For a Fermi liquid, the Korringa relation TiTK"^ = const = {KB? / AT^kB), where R is the 
ratio of nuclear to electronic gyromagnetic ratios, holds but it fails in the presence of 
strong magnetic correlations and/or disorder [53]. In particular, AFM correlations can 
increase x{^. 7^ 0) thus enhancing 1/Ti while not appreciably changing the x(q = 0), 
which determines K [55] . In the present experiment, K was found to vary very littl e with 
X and 1/TiT was found to fit well with a Curie- Weiss form 1/TiT = a + h/ {T + 9). While 
the fitting parameters a, b vary only slightly with x, the "Curie- Weiss temperature" 6 
varies substantially; it is zero at the critical concentration so that the singular part 
of 1/TiT varies as 1/T, as is expected for a second-order antiferromagnetic quantum 
critical point [3J. 

The Kyoto group has combined their results [29j with earlier resistivity [50j and 
structural data to produce the informative phase diagram reproduced here in figure [6] 

The purple x denotes the quantum critical point at Xc = 0.33, which is precisely 
the value of x at which the measured "Curie- Weiss temperature" 9 vanishes. It is also 
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Figure 6. Phase diagram for BaFe2(Asi_a;Pj;)2. The color shading represents 
tlie value of the resistivity exponent in the relation p = po + AT"'. The inverted 
triangles, black dots and squares denote the structural, magnetic (gray region) and 
superconducting (brown region) transitions, respectively. The purple x marks the 
quantum critical point at = 0.33. Courtesy of Y. Matsuda. 

the point where the boundaries of the quantum critical cone, which mark the crossover 
away from a hnear temperature dependence of the resistivity, converge at T = 0. The 
fit in figure |5] shows that it is also the concentration where a critical divergence of the 
effective mass occurs, as determined from thermodynamic measurements. 

5.2. 1111 compounds 

These are the iron pnictides of composition RFePnO, where R is a rare earth, either non- 
magnetic (La) or magnetic (Ce) and Pn is a pnictogen. For LaFeAsO, as the temperature 
is lowered, there is first, at T,, tetragonal to orthorhombic structural second order 
transition followed, within 20 degrees, at T/v, by a second order magnetic transition to 
a collinear antiferromagnetic (AFM) state. In LaFePO, these transitions do not occur: 
there is no magnetic order down to T = and there are conflicting reports [SB] about 
possible superconductivity [57] below 7 K. Neither superconductivity nor magnetic order 
appears in CeFePO. 

In the early work [58], by the Zhejiang group, on LaFeAsi-j^P^jO, superconductivity 
(Tc ~ 10 K) was found just in the region of the disappearance of magnetism and no 
information could be obtained on the magnetic quantum critical point, since the upper 
critical field for the destruction of superconductivity was rather high {Hc2 ~ 27 T). 

The end members of the CeFeAsi_2;Pj;0 series have the following properties: 
CeFeAsO [59] has the successive second order transitions mentioned above (structural 
then magnetic), with Tg 160 K and a Neel temperature T^y of about 135 K. However, 
CeFePO [60] has neither magnetic order nor superconductivity. In this compound 
(x = 1), the Ce 4/ electrons play an important role at low temperature, resulting in 
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heavy fermion properties with a Kondo temperature of about 10 K and ferromagnetic 
(FM) fluctuations from the 4/s. At x = 0, the Ce 4/ electrons are antiferromagnetically 
ordered below 4 K [61] . This order coexists with that of the Fe 3d electrons [62] . Across 
the phase diagram [61], as a function of x, the Ce 4/ order changes from AFM to FM 
in the neighborhood of the Fe 3d quantum critical point at x ~ 0.4, see below. 
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Figure 7. P-doping dependence of the Fe ordered moments and the structural 
orthorhombicity. The proportionaUty between the two is seen as is the vanishing of 
both at the the quantum critical point Xc ~ 0.4. Adapted with permission from |31) . 
Copyright (2010) by the American Physical Society 



At Oak Ridge, neutron scattering studies on CeFeAsi.^^PajO [21] have shown the 
behavior of both T, and T/v as well as the magnitude of the antiferromagnetically 
ordered moment as functions of P doping for As. The separation of the two 
transition temperatures, if any, beyond x = 0.05 is smaller than the experimental 
accuracy. Tg, and the ordered moment all vanish at x = 0.4. Thus, the behavior 
is similar to the blue lines of figure |4j this compound has a magnetic quantum 
critical point free from the interference of superconductivity. As is the case for all 
observations of antiferromagnetism in the pnictides, for the CeFeAsi_^Pi,0 system the 
antiferromagnetic peaks for < x < 0.4 are all commensurate, which supports our use 
of a local moment picture to describe the magnetism. 

The phase diagram (T vs x) of CeFeAsi-ajP^^O has been elaborated using transport 
and thermodynamic measurements by the Zhejiang University group [61]. As is clear 
from figure [8} there is a Fe 3d magnetic quantum critical point at x ~ 0.4, in agreement 
with the neutron scattering results [31]. The very low temperature critical behavior 
(T < 5 K) at X ~ Xc may be distorted by the onset of order in the Ce 4/ moments. In 
contrast to the situation when superconductivity covers the quantum critical point, here 
the interfering phase belongs to different (Ce 4/) degrees of freedom than the ones (Fe 
3d) with which we are concerned. Thus there can be ways of separating the effects. The 
Zhejiang results also confirm the previously observed [60] Ce 4/ heavy fermion behavior 
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Figure 8. Phase diagram of CeFeAsi_a;Pj;0. The inset shows schematically the 
behavior of the Ce 4/ electrons as a function of doping. The competition between the 
Kondo scale (Tk) and the exchange coupling (RKKY) scale J is depicted by the red 
and blue lines. Adapted with permission from [BT. Copyright (2010) by the American 
Physical Society 



at a; = 1. 

The low temperature behavior of various quantities at the critical concentration 
will be influenced by the interplay between the onset of Ce magnetism at about 5 K 
and the Kondo effect with a Kondo temperature of about 10 K. Nevertheless, it may be 
of some interest to examine the behavior of the resistivity and susceptibility at higher 
temperature and x = Xc- These quantities have been measured in [61]. If one tries to 
describe the resistivity by a power law p{T) = po + AT°', then a ~ 1 only above about 
T = 125 K. At lower T, a < 1, which in part reflects the influence of the Kondo effect 
at low temperatures. The experimental susceptibility appears to be dominated by the 
Ce moments and their magnetic transition. 

5.3. Other isoelectronic substitutions 

The substitution of strontium for barium in the 122 compound is also isoelectronic. 
Since the Sr ion is about 17% smaller than the Ba ion, a reduction of the unit cell 
volume V, similar to that for the P for As doping, should be expected. At the critical 
concentration Xc ~ 0.33 for P-doping, V is about 198 |17]. That same volume is 
achieved in Bai_a;Sra;Fe2As2 at x = 0.53. However, the behavior is completely different 
|47j . The antiferromagnetic transition temperature Tjy increases as the unit cell volume 
decreases, rather than decreasing to zero as it does for P-doping. The reason for the 
different behavior lies in the fact that different forms of chemical pressure affect the 
internal structure of the unit cell differently. In figure |9| the relevant parameters are 
shown as hi and h2. Structural studies |17j show that the Fe-pnictogen distance {hi/2) 
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hardly changes at all under Sr-doping so that all of the volume change comes from a 
reduction in h2 — hi, whereas for P-doping, hi decreases substantially, by about 13%. 
Therefore, it is only in the latter case that the Fe-As layer shrinks significantly and of 
course it is just here that the correlation effects develop. 



t 

Ba layer 




Fe 



Figure 9. Schematic of the Ba(FePn)2 (122) crystal structure. The Fe-Pn 
layer thickness is denoted by hi] h2 is one-half the c-axis lattice constant. The a-axis 
lattice parameter is times the Fe-Fe distance. The pnictogen ions (blue dots) lie 
alternately above and below the plane of the figure 

Another possible isoelectronic substitution is Ru for Fe. In Ba(Fei_a;Ru2;)2As2, the 
phase diagram EH ES] for Ru substitution resembles what is depicted in figure [l] 
In spite of initial uncertainty about the valence state of Ru, which does not usually 
occur as 2+, there is evidence that the addition of Ru does not introduce additional 
d-band crossings at the Fermi level and that the total area of hole pockets continues 
to be essentially the same as that of the electron pockets so that n = Uh = rif.. This 
conclusion can be adduced, for example, from both electronic structure calculation |66] 
and ARPES experiment |67]. However, the fact that the 4d orbitals of Ru are larger 
than the 3(i orbitals of Fe leads to substantial changes in the important parameters that 
determine the band structure and the degree of correlation. Thus, the d — p (Ru-As) 
and d — d overlaps increase causing an increase in bandwidth. Interestingly, n appears 
to be about a factor 2 larger for x = 0.35 compared to a; = [6^. In addition the 
effective onsite interaction among the d electrons decreases - both the direct Coulomb 
repulsion U and the Hund's rule coupling. The latter has been argued to be the main 
source of correlation effects and the proximity to a Mott transition in the iron pnictides 
All this implies that Ru substitution takes the system into the weak correlation 
regime, so that the origin of the observed disappearance of magnetism cannot be simply 
ascribed to a bandwidth increase as in the other cases we have discussed, where the 
relevant electron-electron interactions are unaffected by doping. 
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6. Oxychalcogenides 

The chalcogens (Ch) S, Se combine with transition metal atoms (Tm) Mn, Fe, Co to form 
the layered magnetic oxychalcogenide compounds (LaO)2Tm2Ch20. These compounds 
are structurally related to the iron pnictides: Tm2Ch20 layers are separated by (LaO)2 
layers. The former contains a square lattice of Tm ions and when Tm = Fe, it is similar 
to the FeAs layer of a pnictide. The Fe valence is again 2+. In the present context, what 
is interesting about these compounds and analogous ones in which (LaO)2 is replaced by 
(AF)2 where A = Sr or Ba, is that while structurally similar to the undoped pnictides, 
they are all insulating with frustrated antiferromagnetism on a checkerboard lattice 
setting in with Tjv of order 100 K. In a combined theoretical and experimental study 
with collaborators from Los Alamos, Zhejiang, Buffalo, and Rice [HH], we have shown, 
in agreement with earlier work on (AF)2Fe2Ch20 [69], that (LaO)2Fe2Se20 is a Mott 
insulator; this has been confirmed in subsequent studies on this and related compounds 
[701 EH [Z2]- Therefore, these compounds have the PE/KE ratio larger than the value 
for a Mott transition; they are in the Mott insulating state, to the left of the w range 
shown in figure |4j The reason that correlation effects dominate in the oxychalcogenides 
is that the Fe 3d bands are noticeably narrower than those in e.g. LaFeAsO and FeTe, 
reflecting a slight (a few percent) expansion of the Fe-square lattice. Thus, while PE is 
not expected to change much, the KE, measured by the overall bandwidth, is appreciably 
reduced; this drives the material into the Mott state [68]. Thus there is the prospect 
of accessing a Mott critical point in the oxychalcogenides by increasing KE by, for 
example, chemical or external pressure. More generally, the identification of the Mott 
insulator in materials closely related to the parent iron arsenides and chalcogenides but 
with a slightly-expanded Fe-square lattice supports the description, which involves the 
importance of electron-electron correlation, and which forms the basis of the discussion 
in this review. 

7. Conclusion 

We have summarized the theoretical considerations underlying the prediction for a 
magnetic quantum critical point in the parent iron pinictides, and the proposal to 
access it by P doping for As in these materials. The bad-metal nature of the parent 
iron pnictides and chalcogenides places them in the intermediate-interaction region, thus 
in proximity to a Mott transition. In this regime, the spin spectrum comes primarily 
from electronic states away from the Fermi surfaces; we model this by a Ji — J2 spin 
Hamiltonian. At the same time, the electronic spectral weight in the vicinity of the 
Fermi surfaces is tuned by the strength of interactions, and this tuning gives rise to a 
magnetic quantum critical point. 

Isoelectronic doping such as P for As in the parent iron pnictides provides an 
ideal setting to access this quantum critical point. We have gathered the experimental 
evidence that has since been accumulated for such a quantum critical point. In Ce- 
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1111 materials, P doping for As does not induce superconductivity in the measured 
temperature range (above 2K). Neutron scattering measurements of the commensurate 
antiferromagnetic order parameter, combined with transport measurements of the 
temperature-doping phase diagram, make a compelhng case for the existence of 
a magnetic quantum critical point at a P doping about 0.4. In 122 materials, 
superconductivity is observed in the intermediate range of P dopings. Still, the 
combination of transport, NMR, and quantum oscillation measurements provide a strong 
evidence for a quantum critical point in the temperature-doping phase diagram, at a P 
doping of about 0.33, that is however masked by the superconducting dome. 

In carrier-doped iron pnictides, evidence for quantum criticality has been advanced 
[T5| 173] although the case appears not yet as strong as for the isoelectronically-doped 
parent pnictides, which we have discussed here. Moreover, the carrier-doped materials 
are more disordered than the isoelectronically-doped ones, so that the latter are more 
definitively studied by various probes. 

The existence of a quantum critical point in the overall phase diagram of these 
materials has important consequences. The associated quantum criticality serves as 
a likely source of non-Fermi liquid behavior, such as the anomalous temperature 
dependence of electrical resistivity, that has been observed in several different parts 
of the iron pnictide phase diagrams. It provides a proper characterization of the 
magnetic excitations in the iron pnictides and chalcogenides, and this is important 
for understanding the interplay among magnetism, localization, and superconductivity. 
At the same time, it unveils a new setting to advance our general understanding of 
quantum critical phenomena. 

Finally, the isoelectronically-tuned magnetic quantum criticality represents a 
concrete prediction of the strong-coupling approach that has been verified by subsequent 
experiments; this, together with the demonstration of a Mott localization induced by a 
slight expansion of the Fe-square lattice that has also been summarized here, underscores 
the viability of this approach to the microscopic physics of iron-based superconductors. 
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